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■ Abstract 

In the present work, we address a class of Cahn-Hilliard equations characterized by a nonlinear 
diffusive dynamics and possibly containing an additional sixth order term. This model describes 
the separation properties of oil-water mixtures, when a substance enforcing the mixing of the 
phases (a surfactant) is added. However, the model is also closely connected with other Cahn- 
Hilliard-like equations relevant in different types of applications. We first discuss the existence 
of a weak solution to the sixth-order model in the case when the configuration potential of the 
system has a singular (e.g., logarithmic) character. Then, we study the behavior of the solutions 
in the case when the sixth order term is let tend to 0, proving convergence to solutions of the 
fourth order system in a special case. The fourth order system is then investigated by a direct 
approach and existence of a weak solution is shown under very general conditions by means of a 
fixed point argument. Finally, additional properties of the solutions, like uniqueness and parabolic 
regularization, are discussed, both for the sixth order and for the fourth order model, under more 
restrictive assumptions on the nonlinear diffusion term. 

Key words: Cahn-Hilliard equation, nonlinear diffusion, variational formulation, existence theorem. 
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1 Introduction 

This paper is devoted to the mathematical analysis of the following class of parabolic systems: 

u t -Aw = 0, (1.1) 

w = SA 2 u - a(u)Au - ^M|Vu| 2 + f{u) + eu t , (1.2) 

on (0,T) x O, ft being a bounded smooth subset of K 3 and T > an assigned final time. The 
restriction to the three-dimensional setting is motivated by physical applications. Similar, or better 
results, are expected to hold in space dimensions 1 and 2. The system is coupled with the initial and 
boundary conditions 

u|t=o = u , in fl, (1.3) 
d n u = d n w = Sd n Au = 0, on dfl, for t £ (0, T) (1.4) 
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and represents a variant of the Cahn-Hilliard model for phase separation in binary materials. The 
function / stands for the derivative of a singular potential F of a double obstacle type. Namely, 
F is assumed to be +00 outside a bounded interval (assumed equal to [—1, 1] for simplicity), where 
the extrema correspond to the pure states. A physically significant example is given by the so-called 
Flory-Huggins logarithmic potential 

F(r) = (1 -r)log(l -r) + (1 +r)log(l + r) - ^r 2 , A > 0. (1.5) 

As in this example, we will assume F to be at least A-convex, i.e., convex up to a quadratic pertur- 
bation. In this way, we can also allow for singular potentials having more than two minima in the 
interval [—1,1] (as it happens in the case of the oil- water-surfactant models described below, where 
the third minimum appears in relation to the so-called "microemulsion" phase). 

We assume the coefficients 6, e to be > 0, with the case 6 > giving rise to a sixth order model 
and the case e > related to possible viscosity effects that are likely to appear in several models of 
Cahn-Hilliard type (see, e.g., [31]). The investigation of the limits as 6 or e tend to zero provides 
validation of these models as the approximates of the limit fourth order model. 

The main novelty of system (|l.ll) - (ll.2[) is related to the presence of the nonlinear function a 
in (|1.2|) . which is supposed smooth, bounded, and strongly positive (i.e., everywhere larger than some 
constant a > 0). Mathematically, the latter is an unavoidable assumption as we are mainly interested 
in the behavior of the problem when <5 is let tend to and in the properties of the (fourth order) 
limit system 6 — 0. On the other hand, at least in the physical context of the sixth order model, 
it would also be meaningful to admit a to take negative values, as it may happen in presence of the 
"microemulsion" phase (see [35] HH])- We will not deal with this situation, but we just point out that, 
as long as 6 > is fixed, this should create no additional mathematical difficulties since the nonlinear 
diffusion term is then dominated by the sixth order term. 

From the analytical point of view, as a basic observation we can notice that this class of 
systems has an evident variational structure. Indeed, (formally) testing (jl.ip by w, (jl.2l) by u t , taking 
the difference of the obtained relations, integrating with respect to space variables, using the no-Bux 
conditions (|1.4I) . and performing suitable integrations by parts, one readily gets the a priori bound 

±£s(u) + \\Vw\\ 2 L2{n) + s\\u t \\ 2 LHn) = 0, (1.6) 
which has the form of an energy equality for the energy functional 

£ s (u)= f (W + ^|W| 2 +F(u)), (1.7) 

where the interface (gradient) part contains the nonlinear function a. In other words, the system 
(|1.1[) - (|1.2[) arises as the (7J 1 /-gradient flow problem for the functional E$. While the literature on 
the fourth order Cahn-Hilliard model with logarithmic free energy is very wide (starting from the 
pioneering work [35] up to more recent works like, e.g., [T] [23135], see also the recent review [H] and 
the references therein), it seems that potentials of logarithmic types have never been considered in 
the case of a nonconstant coefficient a. Similarly, the sixth order Cahn-Hilliard type equations, which 
appear as models of various physical phenomena and have recently attracted a notable interest in the 
mathematical literature (see discussion below) , seem not to be so far studied in the case of logarithmic 
potentials. 

The sixth order system (|l.lj) - (ll.2j) arises as a model of dynamics of ternary oil- water-surfactant 
mixtures in which three phases occupying a region Q in R 3 , microemulsion, almost pure oil and 
almost pure water, can coexist in equilibrium. The phenomenological Landau-Ginzburg theory for 
such mixtures has been proposed in a series of papers by Gompper et al. (see, e.g., [25J 12S1 120] and 
other references in [H]). This theory is based on the free energy functional (|1.7[) with constant 6 > 
(in general, however, this coefficient can depend on u, see [46), and with F(u), a(u) approximated, 
respectively, by a sixth and a second order polynomial: 

F(u) = (u+ lf{u 2 + h Q ){u- l) 2 , a(u) = g Q +g 2 u 2 , (1.8) 
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where the constant parameters ho, go, gi are adjusted experimentally, gi > and ho, go are of arbitrary 
sign. In this model, u is the scalar, conserved order parameter representing the local difference 
between oil and water concentrations; u = — 1, u = 1, and u — correspond to oil-rich, water- rich 
and microemulsion phases, respectively, and the parameter ho measures the deviation from oil-water- 
microemulsion coexistence. 

The associated evolution system (jl.ip - (jl.2D has the standard Cahn-Hilliard structure. Equa- 
tion expresses the conservation law 

u t + V-i = (1.9) 

with the mass flux j given by 

j = -MVw. (1.10) 

Here M > is the constant mobility (we set M — 1 for simplicity), and w is the chemical potential 
difference between the oil and water phases. The chemical potential is defined by the constitutive 
equation 

S£s(u) . 

w = — \-eu u (1.11) 

ou 

where S£ $^ is the first variation of the functional £s{u), and the constant e > represents possible 
viscous effects. For energy (|1.7[) equation (|1.11[) yields (|1.2j) . We note also that the boundary conditions 
d n u = Sd n Au = are standardly used in the frame of sixth order Cahn-Hilliard models due to their 
mathematical simplicity. Moreover, they are related to the variational structure of the problem in 
terms of the functional (11.71) . However, other types of boundary conditions for u might be considered 
as well, paying the price of technical complications in the proofs. Concerning, instead, the condition 
d n w = 0, in view of (jl.101) . it simply represents the mass isolation at the boundary of f2. 

The system (|l.ll) - (|1.4p with functions F(u),a(u) in the polynomial form (|1.8p . and with no 
viscous term (e = 0) has been recently studied in [¥T] . It has been proved there that for a sufficiently 
smooth initial datum uq the system admits a unique global solution in the strong sense. 

The sixth order Cahn-Hilliard type equation with the same structure as (|l.ip - (11.2[) . 5 > 0, 
polynomial F(u), and negative constant a, arises also as the so-called phase field crystal (PFC) 
atomistic model of crystal growth, developed by Elder et al., see e.g., [35j Q3H IH], and [25] for the 
overview and up-to date references. It is also worth mentioning a class of sixth order convective 
Cahn-Hilliard type equations with different (nonconservative) structure than (jl.ip - (jl.2D . This type of 
equations arise in particular as a model of the faceting a of a growing crystalline surface, derived by 
Savina et al. [UJ (for a review of other convective 4th and 6th order Cahn-Hilliard models see [33] V 
In this model, contrary to (|1.1[) - (|1.2[> . the order parameter u is not a conserved quantity due to the 
presence of a force-like term related to the deposition rate. Such class of models has been recently 
studied mathematically in one- and two- dimensional cases by Korzec et al. (331 1341 135] . 

Finally, let us note that in the case 5 = 0, a(u) = const > 0, the functional (11.7P represents the 
classical Cahn-Hilliard free energy jTTl [18] . The original Cahn-Hilliard free energy derivation has been 
extended by Lass et al. [33] to account for composition dependence of the gradient energy coefficient 
a(u). For a face-centered cubic crystal the following expressions for a(u) have been derived, depending 
on the level of approximation of the nearest-neighbor interactions: 

a(u) — ao + a\u + a^u 2 , (1.12) 

where oo > 0, ai,(i2 6 8 in the case of four-body interactions, 02 = in the case of three-body inter- 
actions, and ai = a,2 = in the case of pairwise interactions. Numerical experiments in |36j indicate 
that these three different approximations (all reflecting the face-centered cubic crystal symmetry) have 
a substantial effect on the shape of the equilibrium composition profile and the interfacial energy. 

A specific free energy with composition dependent gradient energy coefficient a(u) also arises in 
modelling of phase separation in polymers [23) . This energy, known as the Flory-Huggins-de Gennes 
one, has the form (|1.7p with 6 = 0, F(u) being the logarithmic potential (|1.5I) . and the singular 
coefficient 

(1 - u)(l + uj 
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We mention also that various formulations of phase-field models with gradient energy coefficient 
dependent on the order parameter (and possibly on other fields) appear, e.g., in [2l If 4j. 

Our objective in this paper is threefold. First, we would like to extend the result of 41 both 
to the viscous problem (e > 0) and to the case when the configuration potential is singular (e.g., of 
the form l|1.5p ). While the first extension is almost straighforward, considering constraint (singular) 
terms in fourth order equations f p.2j) . in the specific case) gives rise to regularity problems since it 
is not possible, up to our knowledge, to estimate all the terms of equation (|1.2p in L p -spaces. For 
this reason, the nonlinear term f(u) has to be intended in a weaker form, namely, as a selection of 
a nonlinear, and possibly multivalued, mapping acting from V = H (Q) to V . This involves some 
monotone operator technique that is developped in a specific section of the paper. 

As a second step, we investigate the behavior of the solutions to the sixth order system as the 
parameter 6 is let to tend to 0. In particular, we would like to show that, at least up to subsequences, 
we can obtain in the limit suitably defined solutions to the fourth order system obtained setting 5 = 
in (|1.2I) . Unfortunately, we are able to prove this fact only under additional conditions. The reason 
is that the natural estimate required to control second space derivatives of m, i.e., testing (|1.2[) by 
— Au, is compatible with the nonlinear term in Vm only under additional assumptions on a (e.g., if 
a is concave). This nontrivial fact depends on an integration by parts formula devised by Dal Passo, 
Garcke and Grim in |21) in the frame of the thin-film equation and whose use is necessary to control 
the nonlinear gradient term. It is however likely that the use of more refined integration by parts 
techniques may permit to control the nonlinear gradient term under more general conditions on a. 

Since we are able to take the limit S \ only in special cases, in the subsequent part of the 
paper we address the fourth order problem by using a direct approach. In this way, we can obtain 
existence of a weak solution under general conditions on a (we notice that, however, uniqueness is no 
longer guaranteed for 6 = 0). The proof of existence is based on an "ad hoc" regularization of the 
equations by means of a system of phase-field type. This kind of approach has been proved to be 
effective also in the frame of other types of Cahn-Hilliard equations (see, e.g., [B]). Local existence for 
the regularized system is then shown by means of the Schauder theorem, and, finally, the regularization 
is removed by means of suitable a priori estimates and compactness methods. This procedure involves 
some technicalities since parabolic spaces of Holder type have to be used for the fixed point argument. 
Indeed, the use of Sobolev techniques seems not suitable due to the nonlinearity in the highest order 
term, which prevents from having compactness of the fixed point map with respect to Sobolev norms. 
A further difficulty is related with the necessity of estimating the second order space derivatives of u in 
presence of the nonlinear term in the gradient. This is obtained by introducing a proper transformed 
variable, and rewriting (jl.2|) in terms of it. Proceeding in this way, we can get rid of that nonlinearity, 
but at the same time, we can still exploit the good monotonicity properties of /. We note here that a 
different method based on entropy estimates could also be used to estimate Aw without making the 
change of variable, which seems however a simpler technique. 

Finally, in the last section of the paper, we discuss further property of weak solutions. More 
precisely, we address the problems of uniqueness (only for the 4th order system, since in the case 
8 > it is always guaranteed) and of parabolic time-regularization of solutions (both for the 6th and 
for the 4th order system). We are able to prove such properties only when the energy functional £g 
is A-convex (so that its gradient is monotone up to a linear perturbation). In terms of the coefficient 
a, this corresponds to asking that a is a convex function and, moreover, 1/a is concave (cf. |24j for 
generalizations and further comments regarding this condition). If these conditions fail, then the 
gradient of the energy functional exhibits a nonmonotone structure in terms of the space derivatives 
of the highest order. For this reason, proving an estimate of contractive type (which would be required 
for having uniqueness) appears to be difficult in that case. 

As a final result, we will show that, both in the 6th and in the viscous 4th order case, all 
weak solutions satisfy the energy equality (|1.6p . at least in an integrated form (and not just an 
energy inequality). This property is the starting point for proving existence of the global attractor 
for the dynamical process associated to system (|1.1|) - (|1.2|) . an issue that we intend to investigate in a 
forthcoming paper. Actually, it is not difficult to show that the set of initial data having finite energy 
constitutes a complete metric space (see, e.g., [I3J Lemma 3.8]) which can be used as a phase space 
for the system. Then, by applying the so-called "energy method" (cf., e.g., [351 E]), one can see that 
the energy equality implies precompactness of trajectories for t oo with respect to the metric of 
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the phase space. In turn, this gives existence of the global attractor with respect to the same metric. 
On the other hand, the question whether the energy equality holds in the nonviscous 4th order case 
seems to be more delicate, and, actually, we could not give a positive answer to it. 

It is also worth to notice an important issue concerned with the sharp interface limit of the 
Cahn-Hilliard equation with a nonlinear gradient energy coefficient a(u). To our knowledge this issue 
has not been so far addressed in the literature. Let us mention that using the method of matched 
asymptotic expansions the sharp interface limits of the Cahn-Hilliard equation with constant coefficient 
a have been investigated by Pego [42] and rigorously by Alikakos et al. 0. Such method has been 
also successfully applied to a number of phase field models of phase transition problems, see e.g., [16] , 
|15j . In view of various physical applications described above, it would be of interest to apply the 
matched asymptotic expansions in the case of a nonlinear coefficient a(u) to investigate what kind of 
corrections it may introduce to conditions on the sharp interface. 

The plan of the paper is as follows. In the next Section [2] we will report our notation and 
hypotheses, together with some general tools that will be used in the proofs. Section [3] will contain the 
analysis of the sixth order model. The limit 5 \ will then be analyzed in Section^] Section[S]will be 
devoted to the analysis of the fourth order model. Finally, in Section [5] uniqueness and regularization 
properties of the solutions will be discussed, as well as the validity of the energy equality. 

Acknowledgment. The authors are grateful to Prof. Giuseppe Savare for fruitful discussions about 
the strategy of some proofs. 

2 Notations and technical tools 

Let O be a smooth bounded domain of R 3 of boundary r, T > a given final time, and let Q := 
(0,T) x fi. We let H :— L 2 (il), endowed with the standard scalar product (•, ■) and norm || • ||. For 
s > and p 6 [l,oo], we use the notation W s ' p (Sl) to indicate Sobolev spaces of positive (possibly 
fractional) order. We also set H s (fl) := W s ^{&) and let V := H 1 ^). We note by (•, •) the duality 
between V' and V and by || • |[x the norm in the generic Banach space X. We identify H with H' 
in such a way that H can be seen as a subspace of V or, in other words, (V, H, V') form a Hilbert 
triplet. 

We make the following assumptions on the nonlinear terms in (|l.l|l - (|1.2|) : 



The latter condition in (|2.4[) is just a technical hypotheses which is actually verified in all significant 
cases. We also notice that, due to the choice of a singular potential (mathematically represented 
here by assumptions (I2.3|) - (I2.4I) ). any weak solution u will take its values only in the physical interval 
[—1,1]). For this reason, the behavior of a is also significant only in that interval and we have 
extended it outside [—1,1] just for the purpose of properly constructing the approximating problem 
(see Subsection 13.11 below) . 

Note that our assumptions on a are not in conflict with (|1.8p or (|1.12D since these conditions 
(or, more generally, any condition on the values of a(u) for large u) make sense in the different situation 
of a function / with polynomial growth (which does not constrain u in the interval (—1,1)). It should 
be pointed out, however, that the assumptions (|2.ip - (|2.4j) do not admit the singular Flory-Huggins-de 
Gennes free energy model with a(u) given by (I1.13[) . We expect that the analysis of such a singular 
model could require different techniques. 

In (|2.ip , denotes the space of functions that are continuous and globally bounded together 
with their derivatives up to the second order. Concerning /, (|2.3p states that it can be written in the 
form 



aeC t 2 (i;R), 3a,a>0: a<a(r)<aVr£R; 

3a_,a+ G [a, a] : a(r) = a_ Vr < —2, a(r) = a + Vr > 2; 

/eC 1 ^!,!);!), /(0) = 0, 3A>0: f'(r) > -A Vr e (-1, 1); 



(2.1) 
(2.2) 
(2.3) 



lim f(r )r = lim / f - = 

r|->i M-i|/(r)| 



(2.4) 



f(r) = fo(r)-Xr, 



(2.5) 
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i.e., as the difference between a (dominating) monotone part fo and a linear perturbation. By ([2. 3ft - 
(|2.4p . we can also set, for r G (— f , 1), 

Fo(r) := ^ /«,(«) ds and F(r) := F (r) - ^r 2 , (2.6) 

so that F' = f. Notice that Fq may be bounded in (—1,1) (e.g., this occurs in the case of the 
logarithmic potential (|1.5[1 ). If this is the case, we extend it by continuity to [—1,1]. Then, Fq is set 
to be +oo either outside (—1, 1) (if it is unbounded in (— 1, 1)) or outside [—1, 1] (if it is bounded in 
(— 1, 1))- This standard procedure permits to penalize the non-physical values of the variable u and 
to intend fo as the subdifferential of the (extended) convex function Fq : R — > [0, +oo]. 
That said, we define a number of operators. First, we set 

A:V-*V, (Av, z) := / Vv ■ Vz, for v,z E V. (2.7) 

Ja 

Then, we define 

W := {z G H 2 (Q) : d n z = on T} (2.8) 

and recall that (a suitable restriction of) A can be seen as an unbounded linear operator on H having 
domain W. The space W is endowed with the natural i? 2 -norm. We then introduce 

A:W^H, A(z):=-a(z)Az-?-^\Vz\ 2 . (2.9) 
It is a standard issue to check that, indeed, A takes its values in H. 



2.1 Weak subdifferential operators 

To state the weak formulation of the 6th order system, we need to introduce a proper relaxed form of 
the maximal monotone operator associated to the function fo and acting in the duality between V 
and V (rather than in the scalar product of H). Actually, it is well known (see, e.g., [T2] Ex. 2.1.3, 
p. 21]) that fo can be interpreted as a maximal monotone operator on H by setting, for v, £ G H, 

£ = /(,(») inff ^ t(x) = fo(v(x)) a.e. in O. (2.10) 

If no danger of confusion occurs, the new operator on H will be still noted by the letter fo. Corre- 
spondingly, fo is the _ff-subdifferential of the convex functional 

Jo [0,+oo], Fo(v) := f F (v(x)), (2.11) 

where the integral might possibly be +00 (this happens, e.g., when \v\ > 1 on a set of strictly positive 
Lebesgue measure). 

The weak form of fo can be introduced by setting 

£e/o,t»(«) Fo{z)-F (v) foranyzeK (2.12) 

Actually, this is nothing else than the definition of the subdifferential of (the restriction to V of) Fq 
with respect to the duality pairing between V' and V. In general, fo, w can be a multivalued operator; 
namely, fo iW is a subset of V that may contain more than one element. It is not difficult to prove 
(see, e.g., [5J Prop. 2.5]) that, if v s V and fo(v) £ H, then 

Uo(v)} C f 0tW (v). (2.13) 

Moreover, 

if v S V and £ G /o,tu(«) D if, then £ = / (u) a.e. in O. (2.14) 

In general, the inclusion in (|2.13[) is strict and, for instance, it can happen that fo(v) & H (i- e -i v 
does not belong to the i?-domain of fo), while fo yW (v) is nonempty. Nevertheless, we still have some 
"automatic" gain of regularity for any element of fo,w(v): 
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Proposition 2.1. Let v G V, £ G fo,w(v)- Then, £ can be seen as an element of the space .M(ft) — 
C°(ft)' of the bounded real-valued Borel measures on ft. More precisely, there exists T G At (ft), such 
that 

(£, z ) = / z dT for any zeVd C°(ft). (2.15) 
Jn 

PROOF. Let z G C°(ft) n F with z ^ 0. Then, using definition f|2.12f) . it is easy to see that 

(t, z ) = 2\\z\\ L ~ m U,—^ \<2\\z\\ L ~ (n) ((t,v)+T (—^ )-Mv)) (2-16) 

<2||*|| £ oo {n) (Kf J i;>| + |n|(Jb(-l/2) + Fo(l/2))). (2.17) 

This actually shows that the linear functional z h- > (£, z) defined on C°(ft)nV (that is a dense subspace 
of (7° (ft), recall that ft is smooth) is continuous with respect to the sup-norm. Thus, by the Riesz 
representation theorem, it can be represented over C°(ft) by a measure T G VW(ft). I 

Actually, we can give a general definition, saying that a functional £ G V belongs to the space 
V n .M(ft) provided that £ is continuous with respect to the sup-norm on ft. In this case, we can use 
(|2.15[) and say that the measure T represents £ on vVf(ft). We now recall a result (TTJ Thm. 3] that 
will be exploited in the sequel. 

Theorem 2.2. Let v G V, £ £ fo,w{v)- Then, denoting by £ a + £ s = £ the Lebesgue decomposition 
of £, with £ a (£ s ) standing for the absolute continuous (singular, respectively) part of £, we have 



Ca« G i X (ft), (2.18) 

C«(aO = /o(«(a:)) for a.e. i £ ft, (2.19) 

(£,«)- / ^dx^supj / z d£ s , z £ C°(ft), 2 (ft) G [-1, 1]}. (2.20) 

Jn I Jn J 



Actually, in [TT] a slightly different result is proved, where V is replaced by Hq (ft) and, correspondingly, 
M(fl) is replaced by (i.e., the dual of (ft)). Nevertheless, thanks to the smoothness of ft, one 

can easily realize that the approximation procedure used in the proof of the theorem can be extended 
to cover the present situation. The only difference is given by the fact that the singular part £ s may 
be supported also on the boundary. 

Let us now recall that, given a pair X, Y of Banach spaces, a sequence of (multivalued) 
operators T n : X — > 2 Y is said to G-converge (strongly) to T iff 

V (x, y) G T, 3 (x n ,y n ) G T n such that (x n ,y n ) -> (x, y) strongly in X x Y. (2.21) 

We would like to apply this condition to an approximation of the monotone function /o that we now 
construct. Namely, for <j G (0, 1) (intended to go to in the limit), we would like to have a family 
{fa} of monotone functions such that 

fa G C X (M), f' a G i°°(R), / CT (0) = 0, (2.22) 
/cr — > /o uniformly on compact subsets of (—1, 1). (2.23) 

Moreover, noting 

Fa{r) := / / CT (s) ds, for rel, (2.24) 
Jo 

we ask that 

F a (r) > Xr 2 -c, (2.25) 

for some c > independent of a and for all r G M, cr G (0,1), where A is as in (|2.3[) (note that 
the analogue of the above property holds for F thanks to the first of (|2.4|l ). Moreover, we ask the 
monotonicity condition 

F 5j (r) < Fa 2 (r) if cr 2 < <J\ and for all rel. (2.26) 
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Finally, on account of the last assumption (|2.4[) . wc require that 

Vm>0, 3C m >0: f' a (r)-m\f a (r)\>-C m , Vre[-2,2] (2.27) 

with C m being independent of a. Notice that it is sufficient to ask the above property for r 6 [—2, 2]. 
The details of the construction of a family {fa} fulfilling (|2.22p ~ (|2.27p are standard and hence we leave 
them to the reader. For instance, one can first take Yosida regularizations (see, e.g., [T2J Chap. 2]) 
and then mollify in order to get additional smoothness. 

Thanks to the monotonicity property (|2.26p . we can apply [H Thm. 3.20], which gives that 

f a G-converges to f in H x H, (2.28) 
fa G-converges to f , w in V x V' . (2.29) 

A notable consequence of G-convergence is the following property, whose proof can be obtained by 
slightly modifying 7, Prop. 1.1, p. 42]: 

Lemma 2.3. Let X be an Hilbert space, B a , B be maximal monotone operators in X x X' such that 

Ba G-converges to B in X x X' , (2.30) 

as a \ 0. Let also, for any a > 0, v a G X, £ CT G X' such that £ CT G Ba{v a )- Finally, let us assume 
that, for some v G X , £ G X' , there holds 

v a — > v weakly in X, £ CT — > £ weakly in X 1 , (2-31) 
limsup($,a,Va) x < {Z,v) x ■ (2.32) 

<j\0 

Then, £ G B{v). 

Next, we present an integration by parts formula: 

Lemma 2.4. Let u G W H H 3 (Q), £ G V such that £ G fo, w (u). Then, we have that 

(£,Au)>0. (2.33) 

Proof. Let us first note that the duality above surely makes sense in the assigned regularity 
setting. Actually, we have that Au G V. We then consider the elliptic problem 

U a G V, Ua+A 2 Ua + fa{Ua)=U + A 2 U + £ in V' . (2.34) 

Since f a is monotone and Lipschitz continuous and the above right hand side lies in V , it is not 
difficult to show that the above problem admits a unique solution u a G W n H 3 (Q). 

Moreover, the standard a priori estimates for u a lead to the following convergence relations, 
which hold, for some v G V and ( G V, up to the extraction of (non-relabelled) subsequences (in fact 
uniqueness guarantees them for the whole a \ 0): 

Ua — > v weakly in H 3 (Q) and strongly in W, (2.35) 
A 2 Ua — > A 2 v weakly in V, (2.36) 
fa(ua) — > C weakly in V . (2.37) 

As a byproduct, the limit functions satisfy v + A 2 v + ( = u + A 2 u + £ in V'. Moreover, we deduce 
from 

(fa{Ua),Ua) = (u + A 2 U + ^ - U a - A 2 U a ,Ua) , (2.38) 

whence 

lim (fa(ua),ua) = (u + A 2 u + £ - v - A 2 v,v) = ((,v). (2.39) 

a^O v J x ' 

Then, on account of (|2.35|) . (|2.37[) . (|2.29p and Lemma \2. 31 (cf.. in particular, relation (|2.29p ) applied 
to the sequence {fa{va)}, we readily obtain that £ G fo, w (v)- By uniqueness, v = u and ( = 
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Let us finally verify the required property. Actually, for a > 0, thanks to monotonicity of f a 

we have 

< {U(u a ),Au a ) = {u + A 2 u + C- u a - A 2 u a ,Au a ) . (2.40) 
Taking the supremum limit, we then obtain 

< limsup (it + A 2 u + £ - u a - A 2 u a ,Au a ) = (u + A 2 u + £ — u, Au) -liminf (A 2 u a , Au a ) . (2.41) 

Then, using (|2.35l) and semicontinuity of norms with respect to weak convergence, 

-liminf (A*u a ,Au a ) = -liminf II VAitJI 2 < -\\VAu\\ 2 = -(A 2 u.Au), (2.42) 

whence we finally obtain 

< (u + A 2 u + £ - it - A 2 u, Au) = (£, Au) , (2.43) 

as desired. I 

Next, we recall a further integration by parts formula that extends the classical result [121 Lemma 3.3, 
p. 73] (see, e.g., [131 Lemma 4.1] for a proof): 

Lemma 2.5. Let T > and let J : H — >• [0, +oo] a convex, lower semicontinuous and proper 
functional. Let u G ff^O.T; V') n L 2 (0,T;V), T) G L 2 (0,T;V) and let r/(t) G dj(u(t)) for a.e. t G 
(0, T), where dj is the H-subdifferential of J . Moreover, let us suppose the coercivity property 

3fci>0, k 2 >0 such that J (v) > k^vW 2 - k 2 MveH. (2.44) 

Then, the function t H ► J'(u{t)) is absolutely continuous in [0, T] and 

^:J(u(t)) = (u t (t),ri(t)) for a.e. t G (0, T). (2.45) 
at 

In particular, integrating in time, we have 

t 

(u t (r),r)(r)) dr = J{u{t))-J(u(s)) VMe[0,T]. (2.46) 

We conclude this section by stating an integration by parts formula for the operator A. 
Lemma 2.6. Let a satisfy (|2.1I) and let either 

v G H 1 {Q,T;H)r\L 2 {Q,T;W)r\L co {Q), (2.47) 

or 

v G H\0,T;V') nL°°(0,T;FT) n L 2 (0, T; ff 3 (ft)). (2.48) 

Then, the function 

^/^Wl 2 (2.49) 
is absolutely continuous over [0,T]. Moreover, for all s,t G [0,T] we have that 

V«r)),« t (r)) dr= / ^)|V«(t)| 2 - / ^fH|V W ( S )| 2 , (2.50) 

where, in the case (|2.48[) . the scalar product in the integral on the left hand side has to be replaced 
with the duality (v t (r), A(v(r))) . 
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Proof. We first notice that (|2.49D - (|2.50[) surely hold if v is smoother. Then, we can proceed by 
first regularizing v and then passing to the limit. Namely, we define v a , a.e. in (0,T), as the solution 
of the singular perturbation problem 

v rT + aAv rT =V, for ct£ (0,1). (2.51) 

Then, in the case (|2.47l) . we have 

v a G i/ 1 (0,T;^)nL 2 (0,r;i7 4 (O)), (2.52) 

whereas, if (|2.48D holds, we get 

v a G H 1 (0,T;V)nL oo (0,T;H 4 (n)). (2.53) 

Moreover, proceeding as in [20, Appendix] (cf., in particular, Proposition 6.1 therein) and applying 
the Lebesgue dominated convergence theorem in order to control the dependence on time variable, we 
can easily prove that 

v a -> v strongly in H 1 ^, T; H) n L 2 (0, T; W) and weakly star in L°°(Q) (2.54) 

(the latter condition following from the maximum principle), if (|2.47[) holds, or 

v a ^v strongly in if 1 (Q, T; V') n L 2 (0,T; H 3 (fl)) and weakly star in L°°(0,T;W), (2.55) 

if (|2.48|) is satisfied instead. 

Now, the functions v a , being smooth, surely satisfy the analogue of (|2.50l) : 

f (%(r)),^ ( (r)) dr = f ^M| V *v(t)| 2 - f ^^W( S )| 2 , (2.56) 
Js Jn L Ja L 

for all s,t G [0, T]. Let us prove that we can take the limit er \ 0, considering first the case (|2.47D . 
Then, using (|2.54D and standard compactness results, it is not difficult to check that 

A{v a ) -> A(v), (at least) weakly in L 2 (0, T; H). (2.57) 

In particular, to control the square gradient term in A, we use the Gagliardo-Nirenberg inequality 
(cf- M) 

l|Vz|| L 4 (0) < cn\\z\\% 2 \\z\\^l {n) + \\z\\ Vz G W, (2.58) 
so that, thanks also to (12.11) . 

||a'(w (T )|Vw (T | 2 || L2(0!T . ff) < ||a'K)||L=o(Q)||Vw (T |||4 ( Q ) < c\\v a \\ L ^ {Q) (l + \\Av a \\ L 2^ T . H) ) , (2.59) 

and (|237|) follows. Moreover, by ([233]) and the continuous embedding H l {Q,T;H) DL 2 (0,T;FT) C 
C°([0, T]; V), we also have that 

v a -> v strongly in C°([0, T]; V). (2.60) 

Combining (|2~M)) . (|2"3T)) and ([2~pTJ)) . we can take the limit a \ in (|2"36]) and get back (PT50"]) . 
Then, the absolute continuity property of the functional in (|2.49[) follows from the summability of the 
integrand on the left hand side of (|2.50l) . 

Finally, let us come to the case (I2.48|) . Then, (|2.55p and the Aubin-Lions theorem give directly 
(|2.60p . so that we can pass to the limit in the right hand side of (I2.56|) . To take the limit of the left 
hand side, on account of the first (|2.55p . it is sufficient to prove that 

A(v a ) -> A(v) at least weakly in L 2 (0, T; V). (2.61) 

Since weak convergence surely holds in L 2 (0, T; H), it is then sufficient to prove uniform boundedness 
in L 2 (0, T; V). With this aim, we compute 

v(aK)A^ + ^l|V VCT | 2 ) 

= a'{v a )Vv a Av a + a( VrT )VAv a + - ^ |V^| 2 Vi- CT + a {v rT )D 2 v r7 Vv r7 , (2.62) 

and, using (|2.55[) . (|2.1[) . and standard embedding properties of Sobolev spaces, it is a standard proce- 
dure to verify that the right hand side is uniformly bounded in L 2 (0,T;ff) (and, consequently, so is 
the left). This concludes the proof. I 
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3 The 6th order problem 

We start by introducing the concept of weak solution to the sixth order problem associated with 
system (fO)- (fl~4|) : 

Definition 3.1. Let S > and e > 0. Let us consider the 6th order problem given by the system 

u t + Aw = 0, in V', (3.1) 
w = SA 2 u + A(u) +£_ - Xu + eu t , in V' , (3.2) 
C € /o,™(«) (3.3) 

together with the initial condition 

u\ t= o — uo, a.e. in ft. (3.4) 
A (global in time) weak solution to the 6th order problem (|3. 1 [) - (|3 .4[) is a triplet (u,w,£), with 

u e H\0,T;V') nL°°(0,T; W) n L 2 (0, T; iL 3 ^)), eu G i? 1 (0, T; if), (3.5) 

F(«) £L°° (O.TjL 1 ^)), (3.6) 

^L 2 (0,T;y'), (3.7) 

G L 2 (0,T;F), (3.8) 

satisfying ([57T ]) -([57^ |) a.e. in (0,T) together with ([3^4]) . 

We can then state the main result of this section: 

Theorem 3.2. Let us assume (12.1[) - (|2.4p . Let e > and (5 > 0. Moreover, Jet us suppose that 

u £W, F(uo)eL\n), (ti )n G (-1,1), (3.9) 

where (ito)n JS spatial mean of uq. Then, the sixth order problem admits one and only one weak 
solution. 

The proof of the theorem will be carried out in several steps, presented as separate subsequences. 

Remark 3.3. We observe that the last condition in (13.91) . which is a common assumption when dealing 
with Cahn-Hilliard equations with constraints (cf. [35] for more details), does not simply follow from 
the requirement F(uq) G L 1 (r2). Indeed, F may be bounded over [—1, 1], as it happens, for instance, 
with the logarithmic potential (11.51) . In that case, F(uq) G L 1 (fi) simply means — 1 < uo < 1 almost 
everywhere and, without the last (|3.9p . we could have initial data that coincide almost everywhere 
with either of the pure states ±1. However, solutions that assume (for example) the value +1 in a 
set of strictly positive measure cannot be considered, at least in our regularity setting. Indeed, if 
\{u = 1}\ > 0, then regularity (|3.7[) (which is crucial for passing to the limit in our approximation 
scheme) is broken, because f(r) is unbounded for r +1 and £ is nothing else than a relaxed version 
of/(«). 

3.1 Approximation and local existence 

First of all, we introduce a suitably approximated statement. The monotone function /o is regularized 
by taking a family {/ CT }, a G (0, 1), defined as in Subsection 12.11 Next, we regularize uq by singular 
perturbation, similarly as before (cf. (|2.51[1 ). Namely, we take uo jC r as the solution to the elliptic 
problem 

U 0: cr + crAuo,a = u 0l (3.10) 
and we clearly have, by Hilbert elliptic regularity results, 

u ,„eD{A 2 ) Vo-G(G,l). (3.11) 

Other types of approximations of the initial datum are possible, of course. The choice (|3.10p , beyond 
its simplicity, has the advantage that it preserves the mean value. 
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Approximate problem. For a £ (0, 1), we consider the problem 

u t + Aw = 0, (3.12) 
w = SA 2 u + A{u) + f a (u) - Xu+(e + a)u t , (3.13) 
u|t=o = uo,<r, a.e. in fl. (3-14) 

We shall now show that it admits at least one local in time weak solution. Namely, there holds the 
following 

Lemma 3.4. Let us assume (12.1|) - (|2.4|) . Then, for any a e (0, 1), there exist T € (0,T] (possibly 
depending on a) and a pair (u, w) with 

u e H 1 ^, T ; H) n L°°(Q, T ; W) n L 2 (0, T ; I?(A 2 )), (3.15) 
w e L 2 (0,T ;W), (3.16) 

such that (I3.12p ~ d3.13p hold a.e. in (0,Tq) and the initial condition Q3.14p is satishcd. 

Proof. The theorem will be proved by using the Schauder fixed point theorem. We take 

B R := {v S L 2 (0,T ;W)nL 4 (0,T ;W 1A (n)) : |K'||l-(o,t o;W ) + \\v\\mo,T ;W^(n)) < R}, (3.17) 

for To and R to be chosen below. Then, we take u g Br and consider the problem given by (|3.14p 
and 

u t + Aw = 0, in H, (3.18) 
w = SA 2 u + A(u) + fa(u) - Xu + (e + <r)u t , in H. (3.19) 

Then, as u S Br is fixed, we can notice that 

\\Au)\\1 H o,t :H) < 4\Hh(o,T :W) + MtHo,T ;W^(n))) < Q(R)- (3-20) 

Here and below, Q denotes a computable function, possibly depending on cr, defined for any nonneg- 
ative value of its argument (s) and increasingly monotone in (each of) its argument (s). 

As we substitute into (|3.18p the expression for w given by (|3.19p and apply the inverse operator 
(Id+(£ + er)A) -1 , we obtain a parabolic equation in u which is linear up to the Lipschitz perturbation 
/ CT (u). Hence, owing to the regularity (|3.20p of the forcing term, to the regularity (|3.11l) of the 
initial datum, and to the standard Hilbert theory of linear parabolic equations, there exists a unique 
pair (u, w) solving the problem given by (|3.18p - (|3.19p and the initial condition (I3.14p . Such a pair 
satisfies the regularity properties (|3 . 1 5[) - (|3 . 1 6[) (as it will also be apparent from the forthcoming a 
priori estimates). We then note as /C the map such that JC : u *-> u. To conclude the proof we will 
have to show the following three properties: 

(i) K. takes its values in Br; 

(ii) K is continuous with respect to the L 2 (0,T ; W) and the L 4 (0, T ; W 1 ' 4 ^)) norms; 
(in) K. is a compact map. 

To prove these facts, we perform a couple of a priori estimates. To start, we test (|3.18l) by w and 
(|3.19p by Ut (energy estimate). This gives 



dt 



{^Uu\\ 2 + J (F„(u) ±u 2 )^j +(s + a)\\u t \\ 2 + \\Vw\\ 2 



-(A(u),u t )<^\\u t \\ 2 + 7 L\\A(u)\\ 2 (3-21) 



and, after integration in time, the latter term can be estimated using (|3.20p . Next, we observe that, 
thanks to (|2.25l) . we have 



6 -\\Au\\ 2 + f (F„{u) - \u 2 ) > V \\u\\ 2 w - c, (3.22) 
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for some r\ > 0, c > independent of a and for all u in W . Thus, (|3.21[) provides the bounds 

\\ u \\l°°(0,T ;W) + \\ut\\mo,T ;H) + \\^w\\l^(0,T o ;H) < Q(R, T 0, \\uO,a\\w)- (3-23) 

Next, testing (13.19[) by A 2 u and performing some standard computations (in particular, the terms 
(A(u), A 2 u) and (f a (u), A 2 u) are controlled by using p.20[) . Holder's and Young's inequalities, and 
the Lipschitz continuity of f a ), we obtain the further bound 

\\A 2 u\\ L ^ To . H) < Q(R,T , |K CT |k)- (3.24) 

Hence, estimates (|3.23p and p.24p and a standard application of the Aubin-Lions lemma permit to 
see that the range of K, is relatively compact both in L 2 (0,T ;W) and in L 4 (0, T ; W 1,A (SX)). Thus, 
(in) follows. 

Concerning (i), we can now simply observe that, by (|3.23[) . 

\\u\\l2(o,t 0] w) <TQ /2 \\u\\ L ^ {0tTa . w) < T 1/2 Q(i?,T , \\u ^\\ w ). (3.25) 

whence the right hand side can be made smaller than R if To is chosen small enough. A similar 
estimate works also for the L 4 (0, T ; W 1A (n)) -norm since W C W 1A {VL) continuously. Thus, also (i) 
is proved. 

Finally, to prove condition (ii), we first observe that, if {u n } C Br converges strongly to u in 
L 2 (0,To; W) n L 4 (0, Tq; Vy 1,4 (f2)), then, using proper weak compactness theorems, it is not difficult 
to prove that 

A(u n ) -» A(u) weakly in L 2 (0, T : H). (3.26) 

Consequently, if u n (respectively u) is the solution to p.l8[) - p.l9[) corresponding to u n (respectively 
u), then estimates (|3.23[) - p.24j) hold for the sequence {u n } with a function Q independent of n. Hence, 
standard weak compactness arguments together with the Lipschitz continuity of f a permit to prove 
that 

Un = K{u n ) u = JC(u) strongly in i 2 (0, T ; W) ni 4 (0, T ; W 1A (n)), (3.27) 
i.e., condition (ii). The proof of the lemma is concluded. I 



3.2 A priori estimates 

In this section we will show that the local solutions constructed in the previous section satisfy uniform 
estimates with respect both to the approximation parameter a and to the time To. By standard 
extension methods this will yield a global in time solution (i.e., defined over the whole of (0,T)) in 
the limit. However, to avoid technical complications, we will directly assume that the approximating 
solutions are already defined over (0,T). Of course, to justify this, we will have to take care that 
all the constants appearing in the forthcoming estimates be independent of T . To be precise, in the 
sequel we will note by c > a computable positive constant (whose value can vary on occurrence) 
independent of all approximation parameters (in particular of To and a) and also of the parameters e 
and 5. 

Energy estimate. First, integrating p,12p in space and recalling p.lO[) . we obtain the mass con- 
servation property 

(u(t))n = Ka)n = (uo) n . (3-28) 
Next, we can test p.!2p by w, p,13p by u t and take the difference, arriving at 

%r£ atB {u) + IIVHI 2 + (e + <y)\\u t \\ 2 = 0, (3.29) 
at 

where the "approximate energy" £ a j(u) is defined as 




(3.30) 
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Actually, it is clear that the high regularity of approximate solutions (cf. (|3.15|) - (|3.16|> ) allows the 
integration by parts necessary to write (|3.29l) (at least) almost everywhere in time. Indeed, all single 
terms in f|3 . 13[) lie in L 2 (0,T; H) and the same holds for the test function u t . 
Then, we integrate (|3.29l) in time and notice that, by (|2.25[) . 

>»?(*||«||vr+ WAV) -c Vte(0,T). (3.31) 
Consequently, Q3.29[) provides the bounds 

ll«IU«(o,TiV) + 8 1/2 \\u\\ L <*> (0 ,T;W) + (£ + v) 1/2 \\ut\\mo,T-,H) < c, (3.32) 
||Vw|| L 2 (0:T;ff ) < c, (3.33) 
\\F a (u)\\ L co (0 , T . tL i m <c, (3.34) 

where it is worth stressing once more that the above constants c neither depend explicitly on S nor 
on e. 

Second estimate. We test (|3.13[) by u — un, uq denoting the (constant in time) spatial mean of u. 
Integrating by parts the term A(u), we obtain 



6\\Au\\ 2 + / a(u)\Vu\ 2 + / fa(u)(u-u n ) 
Jn Jn 

< (w + Xu- (e + a)u t ,u-un) - [ ?-^-\Vu\ 2 (u ~ u n ) (3.35) 

Jn 2 

and we have to estimate some terms. First of all, we observe that there exists a constant c, depending 
on the (assigned once uq is fixed) value of uq, but independent of a, such that 

f a (u)(u - un) > ^\\fa(u)\\ L i(n) - c. (3.36) 

To prove this inequality, one basically uses the monotonicity of f a and the fact that / CT (0) = 
(cf. Appendix] or [271 Third a priori estimate] for the details). Next, by (|2.2p . the function 
r H> a!{r){r — uq) is uniformly bounded, whence 

~ Jj^\Vu\ 2 {u-un)<c\\Vu\\ 2 . (3.37) 

Finally, using that (wn + Auq, u — uq) = since wn + Xuq is constant with respect to space variables, 
and applying the Poincare-Wirtinger inequality, 

(w + Ait — (e + <r)u t ,u - uq) — (w - wn + A(u - Un) — (e + v)u t ,u - u n ) 
< c\\Vw\\ \\Vu\\ + c\\Wu\\ 2 + c(e + a)\\u t \\ \\Wu\\ 

<c(||VHI + (e + ^)W+l), (3-38) 

the latter inequality following from estimate (|3.32D . 

Thus, squaring (|3.35L using f|3 .36[1 - (|3 .38(1 . and integrating in time, we arrive after recalling 
(j3"32l . ([3331) at 

\\fa(u)\\ L 2 {0tT . L i m) < c. (3.39) 

Next, integrating (|3.13p with respect to space variables (and, in particular, integrating by parts the 
term A(u)), using (13.391) . and recalling (|3.33p . we obtain (still for c independent of e and S) 

IIHU»(o,T;V) ^ c - (3-40) 

Third estimate. We test p,13p by Au. Using the monotonicity of f a and (|2.ip . it is not difficult to 
arrive at 

^^ t \\Vu\\ 2 + 5\\VAu\\ 2 + =\\Au\\ 2 < (Vw + XVu,Vu)+c\\Vu\\t Hn) . (3.41) 
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Using the continuous embedding H 3 / 4 (fl) c L 4 (fl) (so that, in particular, H 7 / 4 (fl) C W 1A {tt)) 
together with the interpolation inequality 

IM| ff 3/4 CO) < ||«fi 8 (n) |Mgf (n) Vv G H^\n), (3.42) 

and recalling estimate (13.321) . the last term is treated as follows: 

c||V«||i4 (n) <c||u||^ a (n) ||«||^ 2 < ^||VAh|| 2 + c05). (3.43) 

Note that the latter constant c(S) is expected to explode as <5 \ but, on the other hand, is indepen- 
dent of cr. Next, noting that 

(Vw + AVu, Vti) < c(||Vu|| 2 + || Vw\\ 2 ), (3.44) 

from (|3.41l) we readily deduce 

IMU 2 (o,T ; H3(n)) < c(6). (3.45) 
A similar (and even simpler) argument permits to check that it is also 

\\A(u)\\ L * {0iT . H) <c(5). (3.46) 

Thus, using (|3.32p . (I3.40[) . (|3.45l) - (|3.46l) and comparing terms in (|3.13[) . we arrive at 

\\U(V)\\Z.>(P,T;V>) < <fy ( 3 - 47 ) 

3.3 Limit cr \ 

We now use the machinery introduced in Subsection l2.1l to take the limit cr \ in (|3.12p ~ (|3.13p . For 
convenience, we then rename as (u a ,w a ) the solution. Then, recalling estimates (|3.32l) - (|3.34[) . (|3.40[) 
and (|3.45l) - (|3.47[) . and using the Aubin-Lions compactness lemma, we deduce 

u a ^u strongly in C°([0,T];i/ 2 - £ (fi)) nL 2 (0,T;iJ 3 - e (fi)), (3.48) 

u a -+u weakly star in H l {0,T;V') n L°°(0, T; VF) n L 2 (0, T; ff 3 (0)), (3.49) 

(e + o> CT)t -> eu t weakly in L 2 (0, T; H), (3.50) 

w a -> w weakly in £ 2 (0, T; V), (3.51) 

UM -> € weakly in i 2 (0,T; V"), (3.52) 

for suitable limit functions u, w, £, where e > is arbitrarily small. It is readily checked that the above 
relations C (|3.48l) in particular) are strong enough to guarantee that 

A(u„) ^A{u), strongly in L 2 (0,T;H). (3.53) 

This allows us to take the limit cr \ in (|3 . 1 (|3 . 14[) (rewritten for u a , w a ) and get 

u t +Aw = 0, in V', (3.54) 
w = 5A 2 u + A{u) + £- Xu + eut, in V, (3.55) 
u\t=o = uq a.e. in f2. (3.56) 

To identify £, we observe that, thanks to (|3.48l) . p.52j) . and Lemma [2.31 applied with the choices of 
X = V, X' = V', B a = /cr, B = f 0lW , v a = u a , v = u and £ CT = f a {u a ), it follows that 

£ 6 /o,™(«). (3-57) 

Namely, £ is identified with respect to the weak (duality) expression of the function /o. This concludes 
the proof of Theorem 13.21 for what concerns existence. 
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3.4 Uniqueness 

To conclude the proof of Theorem l3.21 it remains to prove uniqueness. To this purpose, we write both 
(|3.1j) and (|3.2p for a couple of solutions (u\ , Wi , £1 ) , (u 2 ,w 2 ,£ 2 ), and take the difference. This gives 

u t + Aw = 0, in V', (3.58) 

w = SA 2 u - a(ui)Au - (a(«i) - a{u 2 ))Au 2 - - ^ (|Vm| 2 - |Vu 2 | 2 ) 

- a ' {Ul) ~ a ' {U2) \Vu 2 f + & _ & _ Au + eu u in V, (3.59) 

where we have set (u,w,^) := (ui,Wi,£i) — (1*2, w 2 , £2)- Then, we test (|3.58[) by A~ 1 u, Q3.59P by u, 
and take the difference. Notice that, indeed, u has zero mean value by f|3 . 28|) . Thus, the operator A -1 
makes sense since A is bijective from Wq to Hq, the subscript indicating the zero- mean condition. 

A straighforward computation involving use of standard embedding properties of Sobolev 
spaces then gives 

( - ( Ul )Ati - (o(«i) - a(u 2 ))Au 2 - ^(|V Ul | 2 - |V W2 | 2 ) - °IM^M | V ^ 2 [ 2 , u) 

< Q(|KIU»(o,t ; w), \\u2\\l°°(o.T;W))\\u\\w\\u\\ (3.60) 

and we notice that the norms inside the function Q are controlled thanks to (13.5[) . Thus, also on 
account of the monotonicity of fa. w , we arrive at 

^(l\Hv> + f H 2 ) + ^ll^ll 2 < c||«|k||«|| + A|M| 2 

<c\\u\\%%\\% 3 <^\\A u f + c(S)\\u\\ 2 v ,, (3.61) 

where, to deduce the last two inequalities, we used the interpolation inequality < ||u|| 2 // 3 ||u||™ 3 
(note that (V, H, V) form a Hilbert triplet, cf., e.g., [T31 Chap. 5]) together with Young's inequality 
and the fact that the function || • \\y + \\A ■ || is an equivalent norm on W. Thus, the thesis of 
Theorem 13.21 follows by applying Gronwall's lemma to p. 611) . 

4 From the 6th order to the 4th order model 

In this section, we analyze the behavior of solutions to the 6th order problem as S tends to 0. To start 
with, we specify the concept of weak solution in the 4th order case: 

Definition 4.1. Let 5 = and e > 0. Let us consider the 4th order problem given by the system 



u t + Aw = 0, in V', (4.1) 

w = A(u) + f(u) + su t , in H, (4.2) 

together with the initial condition (|3.4p . A (global in time) weak solution to the 4th order problem 
(j4Tj) - (|4~2)) . (J3T4J) is a pair (u,w), with 

u G -ff 1 (0, T; V') H L°°(0, T; V) fl L 2 (0, T; W), eu G H l (f),T;H), (4.3) 

F(u) e L oo (0,T;L 1 (f2)) ! (4.4) 

f (u) £L 2 (0,T;tf), (4.5) 

weL 2 (0,T;V), (4.6) 

satisfying (|4~T])-([4~2]) a.e. in (0,T) together with ([331). 
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Theorem 4.2. Let us assume (|2.1j) - (|2.4j) together with 

a is concave on [—1,1]. (4-7) 

Let also e > and let, for all 8 € (0, 1), uo.s be an initial datum satisfying (|3.9p . Moreover, let us 
suppose 

u 0> s^u strongly in V, £s(uq,s) ->• £o(uq), where (u )n e (-1, 1). (4.8) 

Let, for any 5 G (0, 1), (u,s, ws, is) be a weak solution to the 6th order system in the sense of Defini- 
tion [37TJ Then, we have that, up to a (nonrelabelled) subsequence of 8 \ 0, 

u 5 ^u weakly star in H 1 {0,T;V')r\L° c {0,T;V)C\L 2 {0,T;W), (4.9) 

eus^eu weakly in H 1 (0,T;H), (4.10) 

w s ^w weakly in L 2 {0,T;V), (4.11) 

Su s ^0 strongly in L 2 (0,T;H 3 (n)), (4.12) 

it -> /o(«) weakiy in L 2 (0, T; V"), (4.13) 

and (u, w) is a weak solution to the 4th order problem. 



Proof. The first part of the proof consists in repeating the "Energy estimate" and the "Second 
estimate" of the previous section. In fact, we could avoid this procedure since we already noted that 
the constants appearing in those estimates were independent of 8. However, we choose to perform 
once more the estimates working directly on the 6th order problem (rather than on its approximation) 
for various reasons. First, this will show that the estimates do not depend on the chosen regularization 
scheme. Second, the procedure has an independent interest since we will see that the use of "weak" 
subdifferential operators still permits to rely on suitable integration by parts formulas and on mono- 
tonicity methods. Of course, many passages, which were trivial in the "strong" setting, need now a 
precise justification. Finally, in this way we are able to prove, as a byproduct, that any solution to 
the 6th order system satisfies an energy equality (and not just an inequality). Actually, this property 
may be useful for addressing the long-time behavior of the system. 

Energy estimate. As before, we would like to test (|3.1[) by wg, (|3.2|l by us,t, and take the difference. 
To justify this procedure, we start observing that ws <E L 2 (0,T;V) by (I3.8[) . Actually, since (|3.1[) is 
in fact a relation in L 2 (0,T; V), the use of wg as a test function makes sense. The problem, instead, 
arises when working on (|3.2[) and, to justify the estimate, we can just consider the (more difficult) 
case e = 0. 

Then, it is easy to check that the assumptions of Lemma 12.61 are satisfied. In particular, we 
have ([ipg} thanks to ([53]) . Hence, (|2~5U1) gives 

Id f 

(us, u A(u s )) = -— / a{u s )\Vus\ 2 , a.e. in (0,T). (4.14) 
Thus, it remains to show that 



<2 



(us,t,8A 2 u s + i s }= (^\Au 5 \ 2 + F(u s j), a.e. in (0,T). 



(4.15) 



To prove this, we observe that 

8A 2 us + is £ L 2 (0, T; V) V6e(0,l). (4.16) 

Actually, we already noted above that ws, A{ug) lie in L 2 (0,T;V). Since us £ L 2 (0,T;V) by (I3.5[) 
and we assumed e = 0, (|4.16l) simply follows by comparing terms in (|3.2I) . Thus, the duality on the 
left hand side of (|4.15p makes sense. Moreover, as we set 

8 



Js(v):= / (-\Av\ 2 +F(v)), (4.17) 
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then a direct computation permits to check that 

SA 2 u s + & G dJs(u s ) a.e. in (0, T). (4.18) 

Indeed, by definition of .ff-subdifferential, this corresponds to the relation 

(8A 2 u 5 + £g,v-ug) < J 5 (v) - Js(u s ) VveH, (4.19) 

and it is sufficient to check it for v £ V since for v G V\H the right hand side is +oo and consequently 
the relation is trivial. However, for v G V, (|4.19|) follows by definition of the relaxed operator fo. w . 
Thanks to (|4. 18[) , (14.151) is a then a direct consequence of inequality (|2.45l) of Lemma 12.51 

Thus, the above procedure permits to see that (any) weak solution (us, wg, £s) to the 6th order 
problem satisfies the energy equality 

^-£g(u(t)) + \\Vw(t)\\ 2 + e\\u t (t)\\ 2 = (4.20) 
at 

for almost all t G [0,T]. As a consequence, we get back the first two convergence relations in (I4.9[) as 
well as (|4.10p . Moreover, we have 

\\Vws\\mo,T:H) <c. (4.21) 

Second estimate. Next, to get (14. lip and (|4.13[) . we essentially need to repeat the "Second estimate" 
of the previous section. Indeed, we see that us — (us)q is an admissible test function in (I3.1[) . However, 
we now have to obtain an estimate of £s from the duality product 

(£s,u s -(ush). (4.22) 

Actually, if £5 = £g,a + is the Lebesgue decomposition of the measure £a given in Theorem 12.21 
then, noting that for all t G [0,T] we have ug(t) G W C C°(Q), we can write 

(Zs(t),u s (t) - (itf)n) = / &,«,(*)(«*(*) ~ (ush) dx + [ (u s (t) - (ush) dfr, a (t). (4.23) 



Next, we notice that, as a direct consequence of assumption 

3^6(0,1): -l+M<(«o,*)n<l-M. V<5g(0,1), (4.24) 

where u is independent of 8. In other words, the spatial means (uo,s)n are uniformly separated from 
±1. Then, recalling (|2.19j) and proceeding as in (|3.36[) . we have 

&,»(*) M*) - (««)n) dx > -Us,a(t)hHn) ~ c, (4.25) 

where c does not depend on 8. 

On the other hand, let us note as d£,s.s = </>s. s d|£,5. s | the polar decomposition of £g, s , where 
\£s.s\ is the total variation of £s,s (c£, e.g., (3H Chap. 6]). Then, introducing the bounded linear 
functional Sg : C° (Q) — > M given by 

Ss(z) := [_z d&,, (4.26) 
using, e.g., [33J Thm. 6.19], and recalling ([2.20j) . we can estimate the norm of S as follows: 

\ts,.\(fi)=j_d\Z6,.\ = \\Ss\\M<n) 

= sup{ J_zd£ s , s , zeC°(n), «(n)c[-l,l]} 

= (^s,ug)— / £s,a,ug = / ug d£g tS = / u S (j)s <s d|& s |, (4.27) 
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where we also used that us £ C°(0). Comparing terms, it then follows 

Us = <ps,s, |&,s| - a.e. in O. (4.28) 

Then, since is clear that 

us = ±l => ,^'i M f , =±1, (4.29) 
\us - (us)n\ 

coming back to (|4.27[) we deduce 

d|6'.s| = / <Ps,s-r^ — 7-4^7 d IC<5,s| < c / 4>s,s(us ~ (us)q) d\£ Sa \ = c [ (us - (u s )n) d£ S ,s- 
! Jn \us-(us)n\ Jq 7sT 

(4.30) 

Here we used again in an essential way the uniform separation property (|4.24|) . 
Collecting (l4T2"3>(|430")) , we then have 

(&,««- (ush) > ~||&,«(t)lk(n) + »7 I d\is, s \-c, (4.31) 



2 

for some c > 0, rj > independent of S. On the other hand, mimicking (I3.35[l - (13.381) . we obtain 

S\\Au s \\ 2 + a\\Vu 5 \\ 2 + (£s,us - (u 5 )n) < c(\\Vw s \\ + e\\u s , t \\ + l), (4.32) 
whence squaring, integrating in time, and using (|4.10l) and (|4.21j) . we obtain that the function 

t >-)• || II £i(fi) + / d|£«5, s (t)| is bounded in L 2 (0,T), independently of S. (4.33) 
Jn 

Integrating now (|3.2I) in space, we deduce 

/ w s = - / a'{us)\Vus\ 2 + ( &-A(u 5 ) n , (4.34) 
Jn z Jn Jn 

whence 

(4.35) 



< c(\\Vu 5 \\ 2 + ||^,«(t)|Ui C „j + J_ d \^s(t)\ + lY 



Thus, squaring, integrating in time, and recalling (|4.2ip and (|4.33[) . we finally obtain (|4.11D . 

Key estimate. To take the limit 5 > 0, we have to provide a bound on A(us) independent of 5. This 
will be obtained by means of the following integration by parts formula due to Dal Passo, Garcke and 
Grim ([HI Lemma 2.3]): 

Lemma 4.3. Let h £ W 2 - co (M.) and z eW. Then, 

h'{z)\Vz\ 2 Az = -- f h"{z)\Vz\ 4 
3 Jn 

+ 1 J h(z)(\D 2 z\ 2 - \Az\ 2 ) h(z)II(Vz), (4.36) 

where II(-) denotes the second fundamental form ofT. 

We then test (|3.2p by Aus in the duality between V' and V. This gives the relation 

~\\S7us\\ 2 + 6\\VAu s \\ 2 + (A(u 5 ),Au s ) + ((s,Au s ) = (Vw s , Vu s ) + X\\Vu s \\ 2 (4.37) 
and some terms have to be estimated. First, we note that 

(A(u s ), Au s ) = (a(u s )Aus + ^^|V^| 2 , Au s ) . (4.38) 
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Thus, using Lemma POl with the choice of h{-) = a'(-)/2, we obtain 

(A(u s ),Au s )= I a{u s )\Au s \ 2 + I f a(u s )(\D 2 u s \ 2 - \Au s \ 2 ) 
Jn Jn 

' f a"(u 5 )|V M5 | 4 + i [ a(u g )n(Vu g ). (4.39) 



Let us now point out that, being T smooth, we can estimate 



a(us)II(yu$) 



< c\\\7us\\ 2 L 2 (r) <u\\Au & \\ 2 w + Cu\\u s \\ 2 , (4.40) 



for small w > to be chosen below, the last inequality following from the continuity of the trace 
operator (applied to Vu) from H s (fl) into L 2 (T) for s £ (1/2, 1) and the compactness of the embedding 
W C H 1+S (fl) for s in the same range. 

Thus, using the concavity assumption (|4.7[) on a and the fact that \u$\ < 1 almost everywhere 
in (0,T) x fl), we get 

(A{us),Au 5 ) >ti\\Au s \\ 2 -c, (4.41) 

for proper strictly positive constants n and c, both independent of S. Next, we observe that, by (13.57)) 
and Lemma T2.41 we obtain (£5, Aus) > 0. Finally, we have 

- (Vwj, Vug) < c||Viy s ||||Vu s ||, (4.42) 

and the right hand side is readily estimated thanks to (|4.10[) and f|4. . 

Thus, on account of (|4.41j) . integrating (|4.37[) in time, we readily obtain the last of (|4.9p as 
well as (|4.12p . Moreover, since — 1 < us < 1 almost everywhere, we have for free 

ll««l|L-((o,r)xn) < !• (4-43) 
Thus, using the Gagliardo-Nirenberg inequality f|2 .58[) . we have also 

u s -> u weakly in i 4 (0, T\ W 1A {Vt)). (4.44) 

This readily entails 

A(u s )-*A(u) weakly in L 2 (0,T;H). (4.45) 
Thus, a comparison of terms in (|3.2j) gives also 

^ -> £ weakly in L 2 (0, T; V')- (4.46) 



Then, we can take the limit S \ in (|3.ip and get (|4.ip . On the other hand, if we take the limit of 
(|3.2[) . we obtain 

io = A(u) + £- \u + eut (4.47) 

and we have to identify Actually, (|4.46p . the strong convergence us — > u in L 2 (0,T; V) (following 
from (|4.9p and the Aubin-Lions lemma) and Lemma 12.31 permit to show that 

?6/o» a.e. in (0,T). (4.48) 

On the other hand, a comparison argument in (14.471) permits to see that £ € L 2 (0,T; H), whence, 
thanks to (|2.14p . we obtain that £(£) — fo(u(t)) S iif for a.e. i S (0,T). This concludes the proof of 
Theorem [ 



5 Analysis of the fourth order problem 

In this section, we will prove existence of a weak solution to Problem (jl.ip ~ p.4p in the fourth order 
case 5 = by means of a direct approach not relying on the 6th order approximation. This will allow 
us to consider a general function a (without the concavity assumption (|4.7l0 . More precisely, we have 
the following 
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Theorem 5.1. Let assumptions (|2.1 |) -(|2.4 p hold, let e > and let 

u eV, FWei 1 ^), (uo)n g (-1, 1). (5.1) 

Then, there exists at least one weak solution to the 4th order problem, in the sense of Definition 14.11 

The rest of the section is devoted to the proof of the above result, which is divided into several steps. 
Phase- field approximation. For a G (0, 1), we consider the system 

Ut + <Jw t + Aw = 0, (5.2) 
w = A(u) + f a (u) — \u + (e + a)u t . (5.3) 

This will be endowed with the initial conditions 

u\ t =o = u 0<a , «;|t=o=0. (5.4) 

Similarly as before (compare with (|3.10jl ). we have set 

U<) )<T + (jA 2 Uo,a = u (5.5) 

and, by standard elliptic regularity, we have that 

u ,<x e H 5 (n) C C 3+a (U) for a e (0, 1/2), d n u , a = d n Au Q:(7 = 0, on T. (5.6) 
Moreover, of course, uo, a uq in a suitable sense as a \ 0. 

Fixed point argument. We now prove existence of a local solution to the phase-field approximation 
by a further Schauder fixed point argument. Namely, we introduce the system 



Ut + crWf + Aw = 0, (5-7) 
w = -a(u)Au - ^7-^|Vu| 2 + f a (u) - Xu + (e + <r)u t , d n u = on T, (5.8) 



which we still endow with the condition (|5.4[) . Here, f a is chosen as in f|2 .22|) . 

Next, we set 

U := {u e C°> 1+a ([0,T o ] x H) : u| t=0 = uo.a, ll«llc'"+= < 2i?} , (5.9) 

where i? := max{l, ||wo,a-|| c i+ot/Qs} and To will be chosen at the end of the argument. It is clear 
that R depends in fact on a (so that the same will happen for To). This dependence is however 
not emphasized here. For the definition of the parabolic Holder spaces used in this proof we refer the 
reader to [371 Chap. 5], whose notation is adopted. Moreover, in the sequel, in place of C°' a ([0, T ] xf2) 
(and similar spaces) we will just write C°' a , for brevity. We then also define 

W := {w e C°' a : w\ t=0 - 0, \\w\\c°.* < R} , (5.10) 

where R is, for simplicity, the same number as in (|5.9[) . 

Then, choosing (u,w) in 14 x W and inserting it in (|5.8[) . we observe that, by the Lipschitz 
regularity of a (cf. (|2.ip ) and standard multiplication properties of Holder spaces, there exists a 
computable monotone function Q, also depending on a, but independent of the time To, such that 

||a(n)|| C o, Q + \\a'(u)\Vu\ 2 \\ c(Ka + \\U(u)\\c^ < Q(R)- (5.11) 

Thanks to [3T, Thm. 5.1.21], then there exists one and only one solution u to (15.81) with the first initial 
condition (|5.4j) . This solution satisfies 

||u||cw« < Q{R)- (5-12) 

Substituting then ut in (|5.7p and applying the same theorem of [3 7) to this equation with the second 
initial condition (|5.4[) . we then obtain one and only one solution w, with 

Mew- < Q{R). (5-13) 
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We then note as T the map such that T : (u, w) H> (u, w). As before, we need to show that: 

(i) T takes its values inWxW; 

(ii) T is continuous with respect to the C° :1+a x C 0,a norm of U x W; 
(in) T is a compact map. 

First of all let us prove (i). We just refer to the component u, the argument for w being analogous 
and in fact simpler. We start observing that, if u € Ui(T(U x W)), III denoting the projection on 
the first component, then 



\\u(t)\\ cam < \\u \\ cam + / \\u t (s)\\ cam ds <R + T Q(R), Vi G [0,T ], (5.14) 

«/ o 

which is smaller than 2i? if To is chosen suitably. 

Next, using the continuous embedding (cf. [37J Lemma 5.1.1]) 

c i,2+a c C 1 / 2 ([0,T ];C 1+a (fi))nC Q/2 ([0,ro];C 2 (fi)), (5.15) 

we obtain that, analogously, 

\\Vu(t)\\c<*m < l|V«o|| -(n) +7 1 o 1/2 ||"|| c1 / 2([0 ,t ];C 1+ °(o)) < iZ + T 1/a Q(i2). (5.16) 

Hence, passing to the supremum for t € [0,To], we see that the norm of u in C , °> 1 + Q can be made 
smaller than 2R if To is small enough. Thus, (i) is proved. 

Let us now come to (Hi). As before, we just deal with the component u. Namely, on account 
of (|5.12p . we have to show that the space C 1,2+a is compactly embedded into C a,1+a . Actually, by 
(|5.15p and using standard compact inclusion properties of Holder spaces, this relation is proved easily. 
Hence, we have (in). 

Finally, we have to prove (ii). This property is however straighforward. Actually, taking 
(u n ,w n ) — > (u,w) in U x W, we have that the corresponding solutions (u n ,w n ) = T(u n ,w n ) are 
bounded in the sense of (|5.12[) - (I5.13[) uniformly in n. Consequently, a standard weak compactness 
argument, together with the uniqueness property for the initial value problems associated to (|5.7[) 
and to (|5.8p . permit to see that the whole sequence (u n ,w n ) converges to a unique limit point (u, w) 
solving (|5.7I) - (I5.8[) with respect to the limit data (u, w). Moreover, by the compactness property 
proved in (in), this convergence holds with respect to the original topology of U x W. This proves 
that (u, w) — T(u,w), i.e., (ii) holds. 

A priori estimates. For any a > 0, we have obtained a local (i.e., with a final time To depending 
on <t) solution to (|5.2p - (|5.3p with the initial conditions (|5.4p . To emphasize the er-dependence, we will 
note it by (u CT , w a ) in the sequel. To let a \ 0, we now devise some of a priori estimates uniform both 
with respect to a and with respect to To. As before, this will give a global solution in the limit and, to 
avoid technicalities, we can directly work on the time interval [0,T]. Notice that the high regularity 
of (u ai w a ) gives sense to all the calculations performed below (in particular, to all the integrations 
by parts). That said, we repeat the "Energy estimate", exactly as in the previous sections. This now 
gives 

||««r|U°°(0,T;V) + \\Fa(Ua)\\L°°(0,T;Li(n)) < C, (5.17) 

(<T + e)V 2 \\u ir , t \\ L , {0tT . B) <c, (5.18) 

CT 1/2 ||Wa|U»(0,T ; H) + \\Vw v \\vno tT .H) < C. (5.19) 

Next, working as in the "Second estimate" of Subsection 13.21 we obtain the analogue of (|3.39p and 
(13~40)1 . 

To estimate fa(u a ) in H, we now test (|5.3p by f^Ua), to get 



£ + Q- d f r u ^ i f ( ^ f , tn N , a'(u a ) 
2 dt 



2 , II J- /„. M|2 



FM+ / (flKl/iK) + < ^Y 1 fa(Uc))\Vu a \ 2 + \\f a {u 



(wv + XuvJ^u,)), (5.20) 
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and it is a standard matter to estimate the right hand side by using the last term on the left hand 
side, Holder's and Young's inequalities, and properties f|5 . 1 T[) and (|3.40p . Now, we notice that, thanks 
to (220 , 

<r)f'M + °^-Ur) > afUr) ~ c\f*(r)\ > §J» - c Vr e [-2, 2], (5.21) 

with the last c being independent of a. On the other hand, for r ^ [—2, 2] we have that a'(r) = by 
(|2.2j) . Hence, also thanks to (I5.17[) . the second term on the left hand side of (|5.20l) can be controlled. 
We then arrive at 

\\UM\\l*(o,T;H) < c. (5.22) 

The key point is represented by the next estimate, which is used to control the second space derivatives 
of u. To do this, we have to operate a change of variable, first. Namely, we set 

cf>(s) := f a 1/2 (r) dr, z a := cf>(u a ) (5.23) 
Jo 

and notice that, by (|2.1[) - (l2.2j) . tfi is monotone and Lipschitz together with its inverse. Then, by (|5.17l) . 

\\z4l~(o,T;V) < c (5.24) 
and it is straighforward to realize that (|5.3p can be rewritten as 



w a = -4>'(u a )Az a + f a o 4> 1 (z a ) - \u<r + (e : + a-)u a ,u d n u a = on T. (5.25) 

By the Holder continuity of u a up to its second space derivatives and the Lipschitz continuity of a and 
a' (cf. (|2.1[) - (|2.2D 'l. — Az CT is also Holder continuous in space. Thus, we can use it as a test function in 
(|5.25p . Using the monotonicity of f„ and </> _1 , and recalling (|5.22[) . we then easily obtain 

IMU 2 (o,tw) < c (5.26) 

Passage to the limit. As a consequence of ([5~T7| - (l5T9l) . (|3~39)) - ([3T40|) and ([5~22l) . we have 

u a ^u weakly star in ^((XT;^') D L°° (0,T;V), (5.27) 

(cr + E)u att -> eu t weakly in L 2 (0, T; H), (5.28) 

faM^f weakly in L 2 {0,T;H), (5.29) 

w a w weakly in L 2 (0,T; V), (5.30) 

u a .t + aw a .t u t weakly in L 2 (0,T;F'): ( 5 - 31 ) 

for suitable limit functions u, w and /. Here and below, all convergence relations have to be intended 
to hold up to (nonrelabelled) subsequences of a \ 0. Now, by the Aubin-Lions lemma, we have 

u a ^u strongly in C°([0,T];H) and a.e. in Q. (5.32) 

Then, (|5.29p and a standard monotonicity argument (cf. Prop. 1.1]) imply that / = f(u) a.e. in 
Q. Furthermore, by (I2.ip - (|2.2p and the generalized Lebesgue's theorem, we have 

a(u a ) — > a(u), a'(u a ) — > a'(u), strongly in L q (Q) for all q E [l,+oo). (5.33) 

Analogously, recalling (|5.24p . z a — 4>(u a ) — > (f>(u) =: z, strongly in L q (Q) for all q £ [1,6). Actually, 
the latter relation holds also weakly in L 2 (0, T; W) thanks to the bound (|5.26l) . Moreover, by (|5.24l) . 
(|5.26l) and interpolation, we obtain 

||Vz CT || L1 o/ 3(Q) <c, (5.34) 

whence, clearly, it is also 

||Vu (T || i io/3 (Q) < c. (5.35) 

As a consequence, being 
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we also have that 

Au a -> Au weakly in L 5/3 (Q). (5.37) 

Combining this with (|5.27[) and using the generalized Aubin-Lions lemma (cf., e.g., [H]), we then 
arrive at 

u a ^u strongly in i 5 / 3 (0,T;W r2 - e ' 5/3 (fi))nC°([0,T];ff 1 ~ £ (^)), Ve > 0, (5.38) 

whence, by standard interpolation and embedding properties of Sobolev spaces, we obtain 

Vtt CT — > Vu strongly in L q (Q) for some q > 2. (5.39) 

Consequently, recalling (|5.33|) . 

a'{u a )\Vu a \ 2 -> a'(u)|Vu| 2 , say, weakly in L^Q). (5.40) 

This is sufficient to take the limit a \ in (I5.3[) and get back (|4.2p . To conclude the proof, it only 
remains to show the regularity (|4.3I) for what concerns the second space derivatives of u. Actually, by 
(|5.26j) and the Gagliardo-Nirenberg inequality (|2.58p , 

z e L 2 (0, T; n L°°(Q) c L 4 (0, T; W 1A (n)). (5.41) 

Thus, we have also u £ L 4 (0,T; W 1,4 (0)) and, consequently, a comparison of terms in (|4.2p permits 
to see that Au £ L 2 (0,T; H), whence (I4.3[) follows from elliptic regularity. The proof of Theorem 15. II 
is concluded. 



6 Further properties of weak solutions 
6.1 Uniqueness for the 4th order problem 

We will now prove that, if the interfacial (i.e., gradient) part of the free energy £$ satisfies a convexity 
condition (in the viscous case e > 0) or, respectively, a strict convexity condition (in the non- viscous 
case e = 0), then the solution is unique also in the 4th order case. Actually, the stronger assumption 
(corresponding to n > in the statement below) required in the non-viscous case is needed for the 
purpose of controlling the nonmonotone part of /(it), while in the viscous case we can use the term 
eut for that aim. 

It is worth noting that, also from a merely thermodynamical point of view, the convexity 
condition is a rather natural requirement. Indeed, it corresponds to asking the second differential of 
£s to be positive definite, to ensure that the stationary solutions are dynamically stable (cf., e.g., [35] 
for more details). 

Theorem 6.1. Let the assumptions of Theorem 15.11 hold and assume that, in addition, 

o"(r)>0, Q"(r)<-«. Vre[-l,l], (6.1) 
where k > if e = and k > if e > 0. Then, the 4th order problem admits a unique weak solution. 
PROOF. Let us denote by J the gradient part of the energy, i.e., 

J :V ->[0,+oo), J{u) := f ^\Vu\ 2 . (6.2) 

Jn 2 

Then, we clearly have 

(J'(u),v)= f (a(u)Vu-Vv+^-\Vu\ 2 v). (6.3) 
Jn \ 2 / 

and we can correspondingly compute the second derivative of J as 

f / q!' (u) I Vzt 1 2 vz \ 
(J"(u)v,z} = J y — ^2 — \-a'(u)vVu- Vz + a'(u)z\7u ■ Vv + a(u)Vv ■ VzJ . (6.4) 
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To be more precise, we have that J'(u) E V and J"(u) E C(V,V) at least for u £ W (this may 
instead not be true if we only have u E V, due to the quadratic terms in the gradient). This is 
however the case for the 4th order system since for any weak solution we have that u(t) E W at least 
for a.e. t G (Q,T). 

From (|6.4p , we then have in particular 

, tii, \ > f ( a" '(u)\Vu\ 2 v 2 „ ,. . „ „ , 9 \ 
(J"(u)v,v) = j [ 2 + 2a'(u)uVu- Vw + a(ii)|Vu| 2 J 

(°(">-^f) |v " 12 - (6 - 5) 

whence the functional J is convex, at least when restricted to functions u such that 

ueW, «(n)c[-l,l], (6.6) 

provided that a satisfies 

a{r)a"{r) - 2a'(r) 2 > VrE [-1,1]. (6.7) 

Noting that 

1\" _ 2{a') 2 -aa" 



(:) 



we have that J is (strictly) convex if 1/a is (strictly) concave, i.e., (|6.1[) holds (cf. also (24l Sec. 3] for 
related results). Note that, in deducing the last inequality in (|6.5p . we worked as if it was a" > 0. 
However, if a"(r) = for some r, then also a'(r) has to be due to (|6.7I) . So, this means that in the 
set {u = r} the first two summands in the right hand side of the first line of (|6.5[) identically vanish. 

That said, let us write both (|3.1[) and p.2[) for a couple of solutions (ui,wi), (^2,^2), and 
take the difference. Setting (it, w) :— (ui,Wi) — (^2,^2), we obtain 

Ut + Aw = 0, (6.9) 
tu = J'( Ul ) - J'(u 2 ) + f{ Ul ) - f{u 2 ) + sut. (6.10) 

Then, we can test (|6.9[) by A~ x u, (|6.10l) by u, and take the difference. Indeed, u = u\ — ui has zero 
mean value by (|3.28p . We obtain 



H 

2dt 



(||u||fr, + e\\u\\ 2 ) + (J'( Ml ) - J'(n 2 ), u> + (/(m) - /(«a),«) - (6.11) 



and, using the convexity of J coming from (16.11) and the A-monotonicity of / (see (|2.3[) ). we have, for 
some function £ belonging to W a.e. in time and taking its values in [—1,1], 



~(\\u\\ 2 v , +s\\u\\ 2 ) +k\\Vu\\ 2 < ~(||«||^ + (J"(0u,u) < X\\u\\ 2 . (6.12) 



Thus, in the case e > (where it may be k = 0), we can just use GronwalPs Lemma. Instead, if e = 
(so that we assumed « > 0), by the Poincare-Wirtinger inequality we have 

A|M| 2 <f||V,i|| 2 + c|M|^, (6.13) 

and the thesis follows again by applying Gronwall's lemma to (|6.12[) . I 



6.2 Additional regularity 

We prove here parabolic regularization properties of the solutions to the 4th order system holding in 
the case of a convex energy functional. An analogous result would hold also for the 6th order system 
under general conditions on a since the bilaplacean in that case dominates the lower order terms (we 
omit the details). 
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Theorem 6.2. Let the assumptions of Theorem 16.11 hold. Then, the solution satisfies the additional 
regularity property 

||u|U~(r,T;H0 + IMU-(r,T;WM(n)) < Qij' 1 ) Vr > 0, (6.14) 

where Q is a computable monotone function whose expression depends on the data of the problem 
and, in particular, on uo. 

Proof. The proof is based on a further a priori estimate, which has unfortunately a formal charac- 
ter in the present regularity setting. To justify it, one should proceed by regularization. For instance, 
a natural choice would be that of refining the fixed point argument leading to existence of a weak 
solution (cf. Sec. [3]) by showing (e.g., using a bootstrap regularity argument) that, at least locally in 
time, the solution lies in higher order Holder spaces. We leave the details to the reader. 

That said, we test (13.11) by w t and subtract the result from the time derivative of (|3.2I) tested 
by lit- We obtain 

~\\Vwf + ~±\\ Ut f + {J"(u)u t ,u t )+ [ /'(«)«? <0. (6.15) 
Then, by convexity of J, 

(j"(u)u t ,U t ) > K\\ut\\ 2 L 2^ T . V y (6.16) 

On the other hand, the A-monotonicity of / gives 

/ /'(«)«* > -A|k||ia ( o,ri*o ( 6 - 17 ) 
Jn 

and, if e = (so that n > 0), we have as before 

- M\u t \\h(Q,T;H) > -g IHIi 2 (0,T;V) ~ c \\ u t \\h(0,T;V) ■ ( 6 - 18 ) 

Thus, recalling the first of (|3.5[) and applying the uniform Gronwall lemma (cf. Lemma 1. 1.1]), it 
is not difficult to infer 

||VH| L o. (TiT .fl ) + £ 1/2 \\ut\\L^(r,T;H) + K\\u t \\ L 2 (TtT . v) < Q{t- v ) Vt > 0. (6.19) 

Next, testing Q3.2p by u — un and proceeding as in the "Second estimate" of Subsection l3.21 but taking 
now the essential supremum as time varies in [r, T] , we arrive at 

IMU~(t,T;V) + ll/(«)IU-(r,Tiii(n)) < Q^ 1 ) Vr > 0. (6.20) 

Thus, thanks to (|6.20p . we can test (I4.2j) by — Az, with z = (j>(u) (cf. (|5.23|) ')- Proceeding similarly 
with Section [5] (but taking now the supremum over [r, T] rather than integrating in time) , we easily 
get (|6.14p . which concludes the proof. I 

6.3 Energy equality 

As noted in Section 21 any weak solution to the 6th order system satisfies the energy equality (|4.20p . 
We will now see that the same property holds also in the viscous 4th order case (i.e., if 6 = and 
e > 0). More precisely, we can prove the 

Proposition 6.3. Let the assumptions of Theorem 15.11 hold and let e > 0. Then, any weak solution 
to the 4th order system satisfies the integrated energy equality 

£b(«(*)) = £b(«o) - / (||v w ( s )|| 2 - £ |M s )|| 2 )d s V*e[o,T]. (6.21) 

Jo 



26 



Proof. As before, we proceed by testing (|4.1[) by w, (|4.2[) by u t and taking the difference. As 
u t G L 2 {0,T;H) and f (u) G L 2 {0,T;H) (cf. (O and (g3])), then the integration by parts 

(/(«),«*) = ^ y a.e. in (0, T) (6.22) 

is straighforward (it follows directly from [T21 Lemma 3.3, p. 73]). 

Moreover, in view of (|4.3I) . assumption (|2.47p of Lemma |2"1)1 is satisfied. Hence, by (|2.50p . we 
deduce that 

f (A(u(s))Ms)) ds= f <M)l\Vu{t)\ 2 - [ ^\Vu Q \ 2 . (6.23) 

Combining (I6.22[) and (|6.23[) . we immediately get the assert. I 

It is worth noting that the energy equality obtained above has a key relevance in the investigation of 
the long-time behavior of the system. In particular, given m G (—1, 1) (the spatial mean of the initial 
datum, which is a conserved quantity due to (|3.28p ). we can define the phase space 

X s<m := {u G V : Su G W, F(u) G L 1 (0), u n = m} (6.24) 

and view the system (both for <5 > and for 5 = 0) as a (generalized) dynamical process in Xg >m . 
Then, (|6.21[) (or its 6th order analogue) stands at the basis of the so-called energy method (cf. [51138]) 
for proving existence of the global attractor with respect to the "strong" topology of the phase space. 
This issue will be analyzed in a forthcoming work. 

Remark 6.4. Whether the equality (|6.21[) still holds in the nonviscous case e = seems to be a 
nontrivial question. The answer would be positive in case one could prove the integration by parts 
formula 

f\u u A(u)+f(u)) = f (^H| Vu (t)| 2 + FKt))) - f (^l| Vu (t)| 2 + FM), (6.25) 
under the conditions 

?iGi/ 1 (0,r ; y , )ni 2 (0,r;Ty)nL oo (g), A(u) + /(«) g l 2 {o,t-v), (6.26) 

which are satisfied by our solution (in particular the latter (|6.26l) follows by a comparison of terms 
in (14.21) . where now e = 0). Actually, if (|6.26l) holds, then both hands sides of (|6 . 25[) make sense. 
However, devising an approximation argument suitable for proving (I6.25|) could be a rather delicate 
problem. 
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